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Abstract 

We present the first results of a two-dimensional MHD simulation of the solar 
magnetogranulation. The medium was assumed to be compressible, gravitationally 
stratified, radiatively coupled, partially ionized, and turbulent. The simulated mag- 
netogranulation evolved over the course of two hours of hydrodynamic (solar) time. 
A surface (magnetic plume-like) mechanism which forms thin magnetic elements was 
Q_[ found to operate during the process of granule fragmentation. The activity of such a 

mechanism suggests that the magnetogranulation can concentrate and intensify the 
global magnetic flux at the boundaries of convective cells and can also form nearly 
vertical compact magnetic flux tubes by involving the weak horizontal photospheric 
field, which may be, in general, of local (turbulent) nature. 

1 Introduction 

When small-scale magnetic structures were discovered in 1973 [5j [23], the investigators 
became aware of the fact that the major part of the magnetic flux on the Sun, outside 
sunspots and pores, is concentrated in these thin elements which are no more than several 
hundred kilometers in size but have magnetic fields of several hundred millitesla. These 
features have been the object of much observational and theoretical research (see the 
reviews by Solanki j2TJ [22], Schiissler [191 120], and Midler [IS]). Of special interest is 
a multidimensional simulation of the interaction between the thermal granulation-scale 
convection and the magnetic field; this simulation can be done by solving a system of 
equations of radiative magnetohydrodynamics for an inhomogeneous, compressible, and 
gravitationally stratified medium. With such time-dependent models the mechanisms of 
formation of compact magnetic elements can be studied in detail, and the adequacy of the 
simulation can be tested as well by comparing the calculated and observed Stokes profile 
parameters of spectral lines. 

The initial magnetic field intensity and configuration are important parameters in 
such calculations. In most studies (e.g., [U [2], EH]) the initial field was assumed to be 
of global nature, and it was therefore homogeneous, vertically oriented, and it filled the 
entire calculation region. The function of the granulation in this case is to concentrate 
the magnetic field near the boundaries of convective cells and intensify the field. The 
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simulation in pfl |2j [IT] lasted only from five to nine minutes of hydro dynamic time, and 
the authors could follow the effect of kinematic and superadibatic mechanisms on the 
formation of magnetic structures only over the lifetime of those convective flows ( "cells" ) 
which existed from the outset. The effect of granule fragmentation (diffusion) on the 
evolution of magnetic structures could not be investigated. 

A different initial field was taken by Brandt and Gadun [3]. It also was homogeneous 
and vertically oriented, but it had the only horizontal scale (70 km), and its initial intensity 
was 8.5 mT. The time of magnetogranulation simulation was quite long (48 min). However, 
the granulation convection was assumed to be quasi- stationary rather than nonstationary 
(two convective flows of the same scale existed in the dimensional simulation region — the 
symmetrical initial conditions permitted them to coexist for an indefinitely long time). It 
was found, nevertheless, that a magnetic configuration with an intensity of about 100 mT 
may form with such a weak local priming field and a strong convective instability. The 
topological pumping and the density gradient effect [21] ensure the initial intensification 
of the magnetic field in deeper layers; the kinematic mechanism concentrates the field 
lines in the vicinity of downflows. The superadiabatic (thermal) effect is the dominant 
factor in the intensification of the magnetic field; it grows when the magnetic field energy 
is close to the equipartition level or exceeds it. In this case the magnetic field suppresses 
the convective flows, and the magnetic configuration situated in the region of downflows 
suddenly begins to descend. 

So, we may state that all published attempts of the direct numerical simulation of 
magnetoconvection ignore a fundamental feature of the solar thermal convection - its 
nonstationarity. There are also some problems in the formulation of initial conditions for 
the magnetic field. 

In this study we present some results of further development of the two-dimensional 
numerical simulation of magnetogranulation. Nonstationary magnetogranulation was sim- 
ulated over the course of two hours of hydrodynamic time with the use of a new concept 
of initial magnetic field. We show that the fragmentation of granules in the presence of 
horizontal photospheric fields results in the formation of new small-scale magnetic struc- 
tures. This mechanism is given the name surface mechanism. We think it important for 
an understanding of magnetogranulation in various objects to establish the role of this 
mechanism in the formation of flux tubes, but here we only describe it, and a detailed 
analysis of the simulation results will be published elsewhere. 

Small-scale magnetic elements are usually called flux tubes j2Tj [22] • Of course, they 
are not tubes in the two-dimensional plane representation, but we arbitrarily reserve this 
traditional name for them. 



2 Two-dimensional radiative magnet ohydro dynamics 



The complete set of equations of radiative MHD for the conditions of compressible gravi- 
tationally stratified turbulent medium had the following form: 
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Figure 1: Initial model and ini- 
tial magnetic field configuration in the 
MHD simulation. Thin lines) isotherms 
for temperatures (from the top down): 
4000, 5000, 6000, 7000, 8000, 9000, 
10000, and 12000 K. Thick lines) mag- 
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Here A is the vector potential, I is a unit tensor, R is the Reynolds turbulence stress 
tensor, g(g) is the free fall acceleration, U = E + B 2 /8ir/p is the total specific energy 
(E = v 2 /2 + e), e is the internal energy, qo is the kinetic energy dissipation of the 
averaged motion at the subgrid level, D m is the magnetic diffusion coefficient, Qr is the 
divergence of the radiation flux vector (radiation heating/cooling); the spatial rectangular 
coordinates x and z describe the simulation region in the horizontal and vertical directions. 
We also assume the magnetic induction B to be equal to the magnetic field intensity H, 
the magnetic permeability being equal to 1. 

We used the ideal gas equation in our calculations; the radiation pressure and possible 
changes in the electron density due to hydrogen ionization and double ionization of 15 
elements were taken into account. The contribution of H~ and H 2 , molecules to the 
ionization equilibrium was taken into account at T < 6000 K. 

The set of equations and the solution of the transfer equation used to determine Qr 
are described in detail in [2 [6j [10] . The radiation effects were treated in our simulation 
in the "grey" approximation. 

Free upper and lower boundary conditions were taken for velocities and thermody- 
namic quantities, i.e., the free inflow and outflow of matter were permitted: the velocity 
components were determined from the condition dV/dz = 0, and the mean internal en- 
ergy and mean density were taken from the initial homogeneous model [61 [10]. Their 
fluctuation profiles at the upper (lower) boundary were the same as in the layer situated 
below (above). The density at the lower boundary was scaled to provide a constant sum 
of the gas, radiation, and magnetic pressures at the horizontal level. The upper and lower 
boundary conditions for the magnetic field were set, assuming that it was of global nature: 
B x = 0, dB/dz = 0. Periodic side conditions were postulated. 

The size of simulation region was 3920 x 1820 km, the spatial step was 35 km, the 
atmospheric layers extended over about 700 km. The internal structure of thin magnetic 
configurations cannot be studied in detail with such a spatial step, but it is sufficiently 
large for the evolution of these configurations to be followed and the calculation costs to 
decrease. 

We calculated two model sequences with a magnetic field and without field, with the 
same initial conditions. The initial model was taken from the time series of HD models 
[7], but the simulation region was diminished. 

Figure [T] shows the initial magnetic field configuration and the initial model. This is a 
bipolar configuration with a field intensity decreasing with height. The mean intensity B 
was 5.4 mT throughout the calculation region. 
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Figure 2: Evolution of the mean magnetic 
field intensity in the simulation region. 



After several similar calculation runs with different initial conditions for the magnetic 
field (including a homogeneous vertically directed field) we concluded that the configu- 
ration shown in Fig. [TJis the most natural one (uniformly distributed). It also provides 
the optimal conditions for the numerical stability of the solution at the initial moment 
of simulation, since the velocity field and the field of thermodynamic quantities, on the 
one hand, and the magnetic field characteristics, on the other hand, turn out not to be 
self-consistent in the initial model. 



Figure [2] displays the average value of B in the simulation region as a function of magne- 
tohydrodynamic (real solar) time. The simulation may be arbitrarily divided into three 
stages - the initial period (to 20 min), which is mainly controlled by the initial conditions 
adopted for the magnetic field, the transition period, when a mutual rearrangement of the 
thermal convection and the magnetic field occurs, and finally, the period of self-consistent 
solutions. 

Figure |3] demonstrates the evolution of the magnetic field and the brightness field 
(visible granulation pattern at a wavelength of 500 nm). The intergranular lanes in Fig. [3^ 
are shown by dark hatching and granules are shown by light hatching. The granules and 
intergranular lanes were distinguished for every individual model with respect to the mean 
intensity level at A 500 nm. The areas with intensities above this level were recognized as 
granules and below this level as intergranular lanes. 

The magnetic field for the continuum formation level (logr^ = 0, tr being the Rosse- 
land optical depth) is shown in Fig. 13b with the field polarities hatched differently. In 
the course of the initial 20 min the concentration and intensification of the magnetic field 
occur in the vicinity of downflows (x = 3600-3800 km). Then a partial reconnection of 
field lines takes place (the reason is, in part, that the initial field was taken as bipolar), 
and local magnetic field concentrations dissipate. Over the period from the 20th to 35th 
minute there are no marked structures in the magnetic field. A clear-cut structuring in 
the field concentration is observed afterwards: the field becomes stronger in intergranular 
lanes. The onset of concentration (formation of new thin magnetic elements) coincides in 
time and space with the moment when granules break up, and the field separates polar- 
ity at that time as well. The intergranular lanes become brighter when the intensity of 
magnetic features exceeds 100 mT. 

The fragmentation of convective flows is due to the action of surface layers, which serve 
as a thermal boundary [18J. For example, the Rosseland absorption coefficient decreases 
by the order of 1000 and the heat capacity by the order of 10 in hot upflows near the 
visible surface in a height interval as small as 100 km [8]. As a result, the hot thermal 
flows are intensely cooled (due mainly to radiation losses) and become unstable. As shown 
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Figure 3: Spatial and temporal evolution of the emerging monochromatic intensity 
I C /<I C > at A 500 nm (a) and the magnetic field intensity B at the logr# = level 
(b) in the 2-D MHD simulation. The intensity I C /<I C > ranges from 0.017 to 1.78 (on 
a linear scale). The field intensity B ranges from -286.7 to +254.0 mT. Dark hatching) 
negative polarity field, light hatching) positive polarity field. 

in [9j [10J E] , the mechanism through which the flows lose their stability is determined 
by their horizontal scale. The flows related to the granules up to l g = 1000-1400 km in 
size dissipate in the main, while larger flows break up. It was found in [8J from 2-D HD 
models that of crucial importance for small-scale flows (l g < 600 km) is the thermal decay 
due to radiation losses predominantly in the horizontal direction. Such flows cool rapidly 
and dissipate. 

Large-scale convective upflows with l g > 1000 km are highly adiabatic, and the major 
cause of their instability is the excess pressure which develops inside them. This pressure 
excess is a fundamental property of the thermal in convection in stratified media (this was 
found even with the simplest polytropic models [1]). Its relative magnitude is especially 
large near the surface (the thermal boundary), and this gives rise to the buoyancy inversion. 
It is the locality where the adiabatically moving matter penetrates into the photosphere 
and is slowed down through the radiation cooling. In addition, large horizontal dimensions 
of convection cells favor the formation of a highly dense passive cloud of matter which 
hangs over the central parts of convection cells. The excessive pressure forces weaken the 
central convective flow and tend to redistribute the convective energy flux in the horizontal 
direction from central regions to periphery. Thereafter the weakened central fragment of 
the large-scale thermal flow cannot provide the balance of forces to support stability near 
the surface, and a stream of matter directed from the surface downwards is formed from 
the dense cloud; the fragmentation thus comes to an end [8]. 

It was shown in [8] that the flows of intermediate horizontal dimensions (l g = 600- 



5 



Time: 1 n 33. m 5 



Magnetic field 




Time: 1"41. m 



500 



-500 





1000 2000 3000 



d 




1000 2000 3000 




1 000 2000 3000 




1 000 2000 3000 




1000 il* 2000 3000 

Horizontal distonce [km] 



1000 2000 3000 

Horizontal distonce [km] 



Figure 4: Left panels) snapshots of velocity field and temperature. Horizontal lines are 
4000, 5000, 6000, 7000, 12000, and 13000 K isotherms (from the top down). Right panels) 
snapshots of magnetic field lines (dotted) and field polarity with dark hatching for negative 
polarity field and light hatching for positive polarity field. 
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Figure 4 (continued): Snapshots of velocity field and temperature as well as magnetic 
field lines and field polarity. Designations are the same. 



1000 km) experience the effect of both mechanisms to the least extent. These theoretical 
results found an observational support [121 H3] ■ 

The fragmentation of the thermal convective flows in the presence of magnetic fields 
maintains the formation of new small-scale magnetic elements. It this case the streams 
of cool matter newly formed near the surface entrap horizontal magnetic field lines, carry 
them downwards, changing the horizontal orientation of the field to the vertical one, and 
intensify the field by extending the field lines. The photospheric horizontal magnetic field 
is quite week, and it meets the freezing-in condition. The field lines are not completely 
reconnected in this case, since the velocity field is not completely symmetrical in the con- 
vective cell which undergoes fragmentation and the cell itself displays clear-cut polarities, 
with the field intensity growing with depth. 

Figure H] illustrate the mechanism described above. The time in this figures corresponds 
to the time scale in Fig. |3l so that the magnetogranulation evolution can be followed 
in parallel in both figures. The first fragment represents the situation obtained in the 
calculation region after l h 33.5 m of modeling. Two flux tubes of opposite polarities are 
concentrated in intergranular lanes. In the right part of the region (x > 2200 km) a well- 
developed convective flow (a granule) is located and suppression of the central fragment 
has already begun in the lower part of this flow. 

A characteristic feature at the moment l h 39 m is the onset of fragmentation in the 
large-scale flow: the surface level between x ~ 3300 km and x ~ 3700 km is sagged, and 
the horizontal magnetic field lines begin to sag as well. 

The fragmentation progresses further from l' l 41 m to l h 43 m : an active downflow de- 
velops due to the convective instability. The field lines are entrained by it. No complete 
reconnection of the entrained field lines is observed, as the convective cell itself has clear- 
cut polarities in the area where the downflow is passing and the field intensity grows with 
depth inside the sell in accordance with the equipartition condition. The field of one 
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polarity only is intensified (the positive polarity in this case), and the field of the other 
polarity weakens (l h 45 m ). 

Finally, in the time interval from l h 45 m to l h 58 m well-developed flux tubes are ob- 
served: one tube (x ~ 3200-3400 km) was formed by the above mechanism, and the 
second tube (x ~ 400-600 km) appeared due to the kinematic effect. However, the su- 
peradiabatic (or thermal) mechanism plays the major role in the field intensification in 
both cases; it is associated with the superconvective instability of downflows in the pres- 
ence of strong magnetic fields. 

The figures demonstrate not only the formation of flux tubes but their dissipation as 
well - the tube found between x = 600 km and x = 800 km at l h 33. 5 m dissipated com- 
pletely by the moment l /i 43 m , and another compact magnetic configuration, between 2000 
and 2200 km, weakens significantly and disappears as a bright feature in an intergranular 
region. 

Compact fields dissipate in our models due to two mechanisms: the internal insta- 
bility of such intensive small-scale features and their weakening caused by reconnection 
of field lines. The internal instability of flux tubes was studied in detail in [13]: it de- 
velops in intensive magnetic structures of a larger scale as a result of potent convective 
instability - this instability gives rise to an intense evaporation of the matter in the sub- 
photospheric layers inside the tube. Under the conditions governing in a gravitationally 
stratified medium this inevitably leads to the expansion of matter along field lines from 
deeper layers and the "reversion" of the superadiabatic effect - downflows are replaced by 
upflows inside the tube. The fields become less compact, and the tube breaks down. 

The other mechanism is clearly recognized in our simulation. It, in essence, is this: 
the tube and the adjacent regions of convective cells may be of different polarities, and 
the kinematic effect, which steadily squeezes the field out of cells, weakens the tube in the 
intergranular region owing to the reconnection of field lines. 

An interesting effect is produced by the common operation of the kinematic mechanism 
and the thermal one in the concentration and intensification of the magnetic field - newly 
formed compact structures with strong magnetic fields and areas with relatively weak 
fields of opposite polarity may be next to one another. This effect can be seen in Figs [3] 
and HI and it is confirmed by observations of magnetic regions with high spatial resolution 

dung. 

Thus, the role of granulation in the formation of small-scale magnetic structures should 
be revised. While it was originally assumed that it only concentrates the global large-scale 
magnetic field in bundles at cell boundaries, our calculations suggest that it is capable of 
forming thin magnetic tubes from the weak photospheric horizontal field. The existence of 
this field is confirmed by observations. In his controversial paper [TJ] Schiissler gives the 
observed distribution of the magnetic field inclination as a function of field intensity. This 
distribution suggests that weak fields are predominantly horizontal in the photospheric 
layers, and so the mechanism described above is quite realistic. 

We obtained a similar distribution from our models by simulating spectral observations 
during the last 25 min of the MHD modeling (Fig. [5]). First we calculated the Stokes 
profiles for the IR Fe I line 1564.8 nm (it is often used in observation programs). Then we 
applied the well-known Stokes diagnostics techniques to the theoretical spectral "scans" 
in order to determine the inclination of the magnetic field to the vertical, 7, and the field 
intensity B. The inclination angle was found from the relation tan 7 = y/Q/V (this is an 
approximate relation, and it is evident from the known expressions for the amplitudes of 
the a-components of the V, Q, and U Stokes profiles), and the field intensity was estimated 
from the distance between the peaks in the red and blue wings of the V profile. In view 
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Figure 5: The magnetic field inclina- 
tion 7 vs. the field intensity B. 

of a lower spatial resolution of the observations, the agreement between the observed and 
theoretical distributions is quite good. This is one more argument for revising the concept 
f the formation of small-scale magnetic elements. 

4 Conclusion 

A direct numerical simulation of the magnetoconvection on the granulation scales suggests 
that the role of granulation in the formation of small-scale magnetic structures should be 
revised. The granulation can concentrate the magnetic field in bundles on the boundaries 
of convective cells and enhance the field intensity in them, but it can also form small-scale 
flux tubes from weak horizontal fields, which may be, in principle, of local nature. 
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